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Abstract 

In a recent paper [Eur. Phys. J C 44, 567 (2005)] we developed a very general formulation 
to take into account explicitly the effects of hydrodynamic flow prohle on the gluonic 
breakup of J/V’’s produced in an equilibrating quark-gluon plasma. Here we apply that 
formulation to the case when the medium is undergoing cylindrically symmetric transverse 
expansion starting from RHIC or LHC initial conditions. Our algebraic and numerical 
estimates demonstrate that the transverse expansion causes enhancement of local gluon 
number density Ug, affects the dependence of the average dissociation rate (P) through 
a partial-wave interference mechanism, and makes the survival probability *S'(pj,) to change 
with Pj, very slowly. Compared to the previous case of longitudinal expansion the new 
graph of ^(py) is pushed up at LHC, but develops a rich structure at RHIC, due to a 
competition between the transverse catch-up time and plasma lifetime. 


PACS numbers: 12.38M 
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1 Introduction 


It is a well-recognized fact that hydrodynamic expansion can signihcantly influence the 
internal dynamics of, and signals coming from, the parton plasma produced in relativistic 
heavy-ion collisions. The scenario oi J/'ip suppression due to gluonic bombardment [l]-[8] 
now becomes very nontrivial because of two reasons: i) the flow causes inhomogeneities with 
respect to the time-space location x and ii) careful Lorentz transformations must be carried 
out among the rest frames of the hreball, medium, and if) meson. In a recent paper [8] this 
nontrivial problem was formally solved by hrst assuming a general flow velocity prohle 
v{x) and thereafter deriving new statistical mechanical expressions for the gluon number 
density ng{x), average dissociation rate (r(a;)), and tjj meson survival probability S{p^) 
at transverse momentum (assuming the meson’s velocity to be along the lateral X 
direction in the hreball frame). 

This general theory was also applied numerically in ref [8] to a plasma undergoing pure 
longitudinal expansion parallel to the collision axis. In such case the kinematics is simple 
because v ■ = 0 and also the cooling is known [9] to occur slowly. When comparison 

was made with the no how situation [7] we found that ng{x) was enhanced, a partial wave 
interference mechanism operated in (r(a;)), and the graph of S{p^) was pushed down/up 
depending on the LHC/RHIC initial conditions. 

The aim of the present paper is to address the following important question: “IT/iaf 
will happen if the general theory of ref [8] is applied to the cylindrically symmetric, pure 
transverse expansion involving tougher kinematics (because v ■ ^ 0) as well as higher 

cooling rate [9]? In Sec.2 below we derive the relevant formulae for statistical observables 
{viz. Ug, (T), S'(py), etc) paying careful attention to the ip meson trajectory and the so called 
catch-up time. Next, Sec.3 presents our detailed numerical work along with interpretations 
concerning (T) and S{p.^). Finally, our main conclusions are summarized in Sec.4. 
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2 Statistical observables 


2.1 Hydrodynamic aspects 


We assume local thermal equilibrium and set-up a cylindrical coordinate system in the 
hreball frame appropriate to central collision. Let x = {r (() z) he Si typical spatial point, 
= it^x) a time-space point, v the fluid 3 velocity, 7 = (1 — the Lorentz factor, 

r the proper time, = ( 7 , 7 ^) the 4 velocity, P the comoving pressure, e the comoving 
energy density, T the temperature, and = {e + P) u^u'^ — Pg^'^ the energy-momentum 
tensor. Then, the expansion of the system is described by the equation for conservation of 
energy and momentum of an ideal fluid 


= 0 , 


( 2 . 1 ) 


in conjunction with the equation of state for a partially equilibrated plasma of massless 
particles 


e = 3P 


a2\g + b2{\ + AJ 


rjnA 


( 2 , 2 ) 


where 02 = &2 = 77r^iVj/40, Nf ^ 2.5 is the number of dynamical quark flavors, A^ 

is the gluon fugacity, and A^ (A^) is the (anti-) quark fugacity. Of course, the gluons (or 
quarks) obey Bose-Einstein (or Fermi-Dirac) statistics having fugacities A^ (or A^). Under 
transverse expansion the fugacities and temperature evolve with the proper time according 
to the master rate equations [ 10 , 11 , 12 ] 


^9tXg + '^^drXg + + ^dri'jT^) 

+ jdrVr + A 

= P 3 (l - A,) - 2 R 2 (^1 - , 

^9tXq + '^^drXq + —dt^-yT^) + 

+ 'ydrVr + A 

V\ V/ 


(2.3) 


(2.4) 
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where the symbols are dehned by 


R 2 — 0.5ng{v(7gg - >qq) , R 3 — 0.5ng{v(7gg - >ggq) (2-5) 

For our phenomenological purposes it will suffice to assume that, at a general instant t in 
the hreball frame, the plasma is contained in a uniformly expanding cylinder of radius 

R = Ri + {t - ti) Ve (2.6) 

where Ri was the radius at the initial instant ti and the expansion speed Ve is a free 
parameter (0 < Ue < !)• In absence of azimuthal rotations the transverse velocity prohle 
of the medium can be parametrized by a linear ansatz 

V = Ve r/R, 0 < r < R (2.7) 

Clearly, |u| vanishes at the origin but becomes at the edge. The (chemical) master 
equations (2.3 - 2.4) are designed to be solved numerically on a computer subject to the 
RHIC/LHC initial conditions stated in Table 1: 

Table 1: Colliding nuclei, collision energy, and initial parameters for the QGP at RHIC(l), 
LHC(l) [13]._ 



Nuclei 

Energy U 

(GeV/nucleon) (fm/c) 

T. 

(GeV) 


A. 

Ri 

(fm) 

RHIC(l) 

Au'^^ 

200 0.7 

0.55 

0.05 

0.008 

6.98 

LHC(l) 

pb208 

5000 0.5 

0.82 

0.124 

0.02 

7.01 


The lifetime or freeze-out time of the plasma is the instant when the temperature 
at the edge falls to T(tijfg) = 0.2 GeV, say. 


2.2 Gluon number density 


For arbitrary flow prohle F, momentum integration [8, eq.(ll)] over a Bose-Einstein distri¬ 
bution function yields the evolving gluon number density 


ng{x) 


16 


TT^ 


iT 


Xg 


oo 

3y-_ 

n=l 


( 2 . 8 ) 
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Since this expression does not depend on the angles of v it has the same structure both for 
the longitudinal and transverse cases. Also, flow enhances the number density compared 
to the no-flow case [7]; e.g. at hxed the enhancement factor 7 becomes 2.3 if I'h'l = 0.9 c. 

2.3 Average ijj dissociation rate 

In the hreball frame (keeping the flow prohle still general) we consider a ip meson of mass 
m^, four momentum p^p, three velocity = pp/Pp’ Lorentz factor 7 ^ = p^/m^. If 
is the plasma 4 velocity measured in the rest frame of ip then we can dehne the useful 
kinematic symbols [ 8 , eq.(30)] 

F = v-v^p, F = - ( 7 p - 1)F 

= 777 (1 --^|fvl) > w = 'j{v-Yv^) 
cos9p,^ = w ■ j {F -Y) /\w\ (2.9) 

where the caps stand for unit vectors. Now, let be the gluon 4 momentum seen in the 
Ip meson rest frame, the cc binding energy, a dimensionless variable, and 

~ the g — Ip breakup cross section according to QCD [14]. Then, 

the mean dissociation rate due to hard thermal gluons [ 8 , eq.(32)] is given by 

X Io{pn) + h{Pn)\v,p\ cos 6 TpW (2.10) 

where we have used the abbreviations 

Cn = ne^w^/T, Dn = ne^\w\/T 
Pn = DnQ^, h{pn) = sinh(p„)/p„ 

h{pn) = cos\i{pn)/pn - sinh(p„)/p 2 ( 2 . 11 ) 

Equation (2.10) demonstrates how (r(a;)) depends on the hydrodynamic flow through 
as well as the angle dpw Retaining only the n = 1 term and picking-up the dominant peak 
contribution from = 10/7 we arrive at the useful approximation 

(f (x)) oc Ag7pi7 

H = \h{D^Ql) + h{D^Ql) \ \ cosOp^J (2.12) 
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in which a partial wave interference mechanism operates due to the anisotropic cos 9^^ 
factor. Numerical consequences of (2.10) relevant to transverse flow will be discussed later 
in Sec.3.1. 

2.4 J/'ip Survival probability 

In this section we shall consider pure transverse flow parametrized by (2.7) and the ip 
meson moving in the lateral X direction with velocity 0 0) appropriate to the 

mid-rapidity region in the hreball frame. Suppressing the coordinate the production 
conhguration of tp meson is called (tj,f^^ = (ti,r^,(p^'^ and the general trajectory after 
time duration A is termed {t,f^) = (t,r^,(p^^'^ such that 

tj ti , A t tj 

= + (2-13) 

where Tp ~ 0.89 fm/c is the proper formation time [15] of the cc bound state. This 
transverse trajectory will hit the edge R = Rj + Ve X of the radially expanding cylinder 
{cf.{2.6)) at the catch-up instant t* after duration A* such that 

I Ri + VeA* |" = I f^ + v^A* f 

so aA*^ + 2f3A* - p = 0 

with a = — vl, 11 = R] — r^p^ 

(3 = rlpVpj cos (plp - R^Ve (2-14) 

If the quadratic in A* has real roots we pick up that which is positive and smaller; but if 
both roots are imaginary then catch-up cannot occur. The time interval of physical interest 
becomes 

ti <t< tii, tn = min {ti + A*, (2.15) 

This formula is quite different from that derived in the case of longitudinal flow [8, eq.(48)]. 
As the time t progresses the dissociation rate (2.10) must be evaluated on the tp meson 
trajectory itself, implying that we have to set at a general instant 

f=f^, v = Ver^/R 

F = v-v^= {rip cos <plp + v^ a) (2.16) 
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in the kinematic relations (2.9). Clearly, the notation (F) of (2.10) becomes equivalent to 




(2.17) 


depending parametrically on the production conhguration 0^. Then, by using the 
radioactive decay law without recombination and averaging over the cross sectional area 
Ai = (at the production instant) we arrive at the desired survival probability 


S{pt) = d\{R] - rlp^)e ^- 4 ^) 

W = f dt f[t], d^rlf^ = dr^ 

Jti 


(2.18) 


Here no information is needed about the length Lj of the cylindrical plasma in contrast 
to the case of longitudinal flow [8, eq.(52)] where the averaging had to be done over the 
volume Vj = nBjLi. 


3 Numerical results 

3.1 Curves of dissociation rate 

The exact formula (2.10) of (T) is a very complicated function of t as well as several 
kinematic parameters dehned jointly by (2.9, 2.11, 2.16) but a feeling for its behaviour 
can be obtained in the extreme nonrelativistic (|n|/c —0) and ultrarelativistic {\v\/c ^ 1) 
limits. For simplicity, suppose at the instant tj a special ip was formed almost at the edge 
Ri of the cylinder with 0^ being the angle between the tjj position vector and velocity 
vector. Then the kinematic relations (2.16, 2.9) yield 

V = F = v^ cos Pip = ±v^ 

w = '^{v-YvpPj ='^e{Fv^-Y)vp, (3.1) 

where the +, — signs correspond to cosplp = ±1, he., to plp = 0,7r, respectively. Thus we 
have the parallel or anti-parallel property 

w II -0^, cos = +1 ii Pip = 0 k. Y < 

w II -vpj, cosppj^ = -1 if 0^ = TT or F > (3.2) 
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Figure 1: The variation of the modihed rate (F) as a function of temperature at different 
transverse momenta for the transverse flow velocity v = 0.2 c for (a) 0^ = 0 and (b) 0^ = tt, 
respectively. 

Figures 1-4 depict the corresponding exact curves of (F) computed from (2.10) based on 
the LHC initial conditions of Table 1. We now proceed to interpret these graphs using the 
approximate estimate (2.12). 

Interpretation; i) At fixed (Py,0^,t^J the steady increase of (F) with T in Figs.l - 2 is 
caused by the growing exp{ —(Ci =F Di)Q^} factors occurring in the estimate (2.12). ii) 
At hxed value of (T, 0^, = 0.2) corresponding to nonrelativistic flow the variation of (F) 

with in Figs.3a,b is more intricate. At 0^ = 0 in Fig.3a there is a broad enhancement 
of (F) for low Pj, < I GeV; this is because, hrstly low speeds of the 0 and plasma can 
compete, and secondly constructive interference occurs between Jq and Ji in the estimate 
(2.12) for cos9^w = +1 (c/.(3.2)). On the other hand, at 0^ = tt in Fig.3b our (F) decreases 
monotonically with p^ throughout; this is due to the fact that, since cosO^pw = — 1 now 
(c/.(3.2)), the interference between Jq and A becomes destructive, iii) At hxed values of 
(T, 0^,u^ = 0.9) corresponding to ultrarelativistic how similar trends with respect to Pj, 
are again explained in Figs.4a, b except for the fact that the steady rise of (F) with in 
Fig.4a is caused mainly by the yp coefficient present in the estimate (2.12). 
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phi=0, v_e=0.9 



phi=pi, v_e=0.9 



Figure 2: The variation of the modihed rate (F) as a function of temperature at different 
transverse momenta for the transverse flow velocity v = 0.9 c for (a) 0^ = 0 and (b) 0^ = tt, 
respectively. 

3.2 Curves of survival probability 

For a chosen creation conhguration of the 0 meson the function W was hrst computed from 
(2.18) and then S{p^) was numerically evaluated. Figures 5a and 5b show the dependence 
of *S'(pj,) on Py corresponding to the LHC and RHIC initial conditions, respectively (for 
two choices of the transverse expansion speed v^). For the sake of direct comparison, we 
also include our earlier results based on no flow [7, eq.(25)] and longitudinal expansion [8, 
eq.(52)] (starting from two possible lengths Lj of the initial cylinder). We now turn to a 
physical discussion of these graphs. 

Interpretation; In every scenario of gluonic dissociation the function W = dt (F) 
depends on via the integrand (F) as well as the limits Three interesting cases 

may now be distinguished: 

No flow case [7]: Here cooling of the plasma is simulated through the master rate equations 
but the existence of the explicit flow prohle is ignored. Then (F) decreases monotonically 
with p^ because of a destructive interference between the Jq and Ji terms. Also, the time- 
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(a) 


phi=pi, v_e=0.2 



(b) 


Figure 3: The variation of the modified rate (f) as a function of transverse momentum for 
different values of temperatures for a transverse flow velocity v = 0.2 c for (a) 0^ = 0 and 
(b) 0^ = TT, respectively. 

span tjj — is shortened as the speed of the tjj meson increases. Consequently, the survival 
probability called 5'o(pj,) grows steadily with as shown by the solid lines in Figs.5a,b. 

Longitudinal expansion case: Here an extra parameter appears namely the length Lj of 
the initial cylinder. For nonrelativistic flow emanating from short length Lj = 0.1 fm the 
(F) values are somewhat reduced compared to the no flow case (due to Jq, h destructive 
interference) though the time span —tj remains unaltered, so that the survival probability 
called 5'||(pj,) is pushed slightly upwards in Figs.5a,b. But for relativistic flow emanating 
from longer length Lj = 1 fm the shifts of the 5'||(py) curve occurs in mutually opposite 
directions at LHC and RHIC (due to the different initial temperatures generated therein). 

Transverse expansion case: Here the extra parameter involved is the transverse expansion 
speed which together with 0^ and T control the p^-dependence of the function W. For 
0^ = 0 the (F) values in Figs.3a, 4a exhibit enhancement/rising trend on the lower p^ side; 
such tfj mesons contribute sizably to W but little to e“^. On the othe hand, all curves of 
(F) in Figs.3 - 4 flatten-off to nearly constant values on the higher p^ side; such ip mesons 
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1 


phi=0, v_e=0.9 



(a) 


phi=pi, v_e= 0.9 



Figure 4: The variation of the modihed rate (F) as a function of transverse momentum at 
different values of temperatures for a transverse flow velocity v = 0.9 c for (a) 0^ = 0 and 
(b) 0^ = TT, respectively. 

contribute substantially to especially for low temperatures. Therefore, the transverse 
survival probability 5'^(Py) becomes nearly Pj,-independent (or very slowly varying) in 
Figs.5a, 5b in sharp contrast to the longitudinal case. For explaining the magnitude of 
the ratio S^{p^)/S\\{p^) we consider the temporal scenario dealing with the limits of the 
integration. 

Temporal scenario: It is known that transverse expansion of a quark-gluon plasma 
produces cooling at a faster rate compared to longitudinal expansion so that the inequality 
fjjfe < holds on the corresponding lifetimes. At LHC the transverse cooling is so 
fast that, for most ip mesons of kinematic interest we have tn = in the dehnition 
(2.15). The time-span tjj —t^ is, therefore, much smaller compared to the longitudinal case 
implying thet S±{p^) > 5'||(pj,) in Fig.5a. Clearly this property at LHC is devoid of any 
rich structure. 

However, at RHIC let us divide the tp meson kinematic region into two parts. For 
slower mesons having < 5 GeV the catch-up time tj + A* in (2.15) exceeds the lifetime 
so that tn = again, i.e., S'^(py) > 5'||(pj,) in Fig.5b for p^ < 5 GeV. Next, for faster 
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LHC(l) 



RHIC(l) 



(b) 


Figure 5: The survival probability of J/'^ in an equilibrating parton plasma at (a) LHC(l) 
and (b) RHIC(l) energies with initial conditions given in Table 1. The solid curve 5'o(pj,) 
is the result of [7], i.e., in the absence of flow while the dotted and dashed curves represent 
the S\\{pj,) when the plasma is undergoing longitudinal expansion with the initial values 
of the length of the cylinder Li = 0.1 fm and 1 fm, respectively [8]. The dot-dashed and 
and double dot-dashed curves depict the S±{p,^) when the system is undergoing transverse 
expansion with the expansion speed = 0.1 and 0.9, respectively. 

mesons having > 5 GeV the reverse inequalities hold making S±{pr^) < Sflpr^) in Fig.Sb. 
Clearly, the rich structure in S±{pj,) at RHIC arises from a mutual competition between 
the catch-up time and lifetime. 


4 Conclusions 

a) In this work we have applied our general formulation [8] of hydrodynamic expansion to 
study the effect of explicit transverse flow prohle on the gluonic break up of J/'^’s created 
in an equilibrating QGP. The formalism in Sec.2 and numerical results of Sec.3 are new 
and original. 

b) Equation (2.8) shows that, at specihed fugacity A^, the effect of transverse flow is to 
increase the gluon number density Ug. This was also the case with longitudinal flow. 
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c) Our expressions (2.10, 2.12) of the mean dissociation rate (F) involves hyperbolic 

functions as well as partial wave interference mechanism (controlled by the anisotropic 
cos6^u] factor). In addition, knowledge about a nontrivial kinematic function F (c/.(2.16)) 
is needed for interpreting the variation of (F) with T, p^, in Figs.l - 4. In contrast, 

for longitudinal flow the treatment of (F) was easier because F = 0 there. 

d) There are several features of contrast between the transverse and longitudinal sur¬ 
vival probabilities denoted by 5'j_(pj,) and 5'||(py), respectively. Due to the geometry of 
production conhguration our S^{p^) contains a double integral (2.18) whereas >S'||(pj,) con¬ 
tains a triple integral. Next, due to the flattening-off trend of (F) with increasing p^ our 
S±{Pj,) becomes roughly ^^.-independent (or slowly varying) in Figs.5a, 5b whereas 5'||(pj,) 
rises rapidly. Finally, the quick cooling rate at LHC makes S±{pj,) > 5'||(pj,) at all p^, of 
interest in Fig.5a whereas a competition between the catch-up time and lifetime generates 
richer structure at RHIC in Fig.5b. 

e) We conclude with the observation that the held oi J/ip suppression due to gluonic 
break up continues to be a research area of great challenge. In a future communication 
we plan to study the effect of asymmetric how prohle arising from noncentral collision of 
heavy ions at hnite impact parameter b. 
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